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Abstract 

This paper is a first report on the cooperative effort in 
helicopter Flight Mechanics Modeling being carried out 
under the agreements of the United States-Israel Memo- 
randum of Understanding. It presents two different mod- 
els for the AH-64 Apache Helicopter which mainly differ 
in their approach to modeling the main rotor. The first 
model, referred to as “BEMAP” (Blade-Element Model 
for the APache), was developed at the Aeroflightdynam- 
ics Directorate, Ames Research Center, and is the only 
model of the Apache to employ a direct blade-element ap- 
proach in calculating the coupled flap-lag motion of the 
blades and the rotor force and moment The second model 
was developed at the Technion-Israel Institute of Tech- 
nology and uses a harmonic approach to analyze the ro- 
tor. This approach allows different levels of approxima- 
tion, ranging from “first-harmonic” (similar to a tip-path- 
plane model) to “complete high-harmonics” (comparable 
to a blade-element approach). The development of the two 
models is outlined and the two are compared using avail- 
able flight test data. 

Introduction 

Mathematical models intended for flight mechanics appli- 
cations are well thought out compromises between sim- 
plicity and accuracy. Generally speaking, the more so- 
phisticated a model is, the more accurately it can match 
the responses of the actual flight vehicle. Sophistication, 
however, brings with it increased costs both in develop- 
ment and in eventual use. For example, models such as 
CAMRAD (Ref 1) are highly sophisticated and multi- 
disciplinary but are cumbersome to use for parametric 
studies in handling qualities and are unsuitable for real- 
time applications. Lack of sophistication, on the other 
hand, can lead to unacceptable inaccuracies. For example, 
stability-derivative-type models, such as TMAN (Ref 2), 
while very simple for real-time applications, can lead to 


erroneous conclusions since they are only applicable to a 
very limited region of the flight envelope, and are not very 
accurate even there. A careful determination of the level 
of sophistication needed to achieve the required accuracy 
is therefore necessary, especially if the model is intended 
for real-time in addition to nonreal-time use. 

The single most important module of any 
component-type helicopter mathematical model is the 
main rotor. This is true not only because the rotor is sin- 
gularly responsible for almost all the forces and moments, 
but also because all other components are significantly af- 
fected by the rotor. The sophistication and accuracy of the 
rotor module, therefore, largely determines the sophistica- 
tion and accuracy of the entire model. There are several 
approaches currently used in flight mechanics for model- 
ing the main-rotor. These include 1) the tip-path-plane ap- 
proach, 2) the rotor-map approach, and 3) the direct blade- 
element approach. All of the above approaches are in- 
herently similar in that they all start with a strip-theory 
modeling of each blade. In the tip-path-plane approach, 
such as that used in ARMCOP (Ref 3), the equations of 
motion are transformed to a non-rotating frame using the 
multiblade-coordinate transformation and solved analyti- 
cally. As a consequence, only very simple linear aerody- 
namics are considered and effects such as compressibility, 
blade stall, and reverse flow are neglected. 

The rotor-map approach, such as used in FLYRT (Ref 4), 
was initially developed in order to allow real-time opera- 
tion of a blade-element rotor. In this approach, a nonreal- 
time blade-element model is run off-line for a great num- 
ber of flight conditions and the results recorded in quasi- 
static look-up tables. The tables are then used by the real- 
time rater module to instantly determine the quasi-static 
rotor forces, moments, and attitudes based on the input 
parameters. Rotor dynamics are then added to the quasi- 
static results to complete the rotor output. This approach is 
also restrictive with regard to modeling secondary effects 
such as compressibility, stall, and reverse flow. 
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Finally, in the direct blade-element approach, such as in 
GEN HEL (Ref 5), the intermediate step of generating 
output tables is eliminated thanks to the power and speed 
of the current computers. This allows easy access to the 
calculations at the elemental level which in turn makes it 
easier to employ sophisticated aerodynamic theories and 
account for the secondary effects in detail. 

Researchers from the U.S. Army Aeroflightdynamics Di- 
rectorate and the Teehnion-Israel Institute of Technology 
have been working cooperatively under the agreements of 
the United States-Israel Memorandum of Understanding 
(MOU) on “Helicopter Flight Control and Display Tech- 
nology” to evaluate alternative methods of formulating ro- 
tor system models for rotorcraft flight-mechanics simula- 
tion. 

While both research groups have used a strip-theory ap- 
proach, the implementations are quite different and each 
has distinct advantages and limitations. The U.S. Army 
has developed a model based on symbolically generated 
exact equations of motion and numerical summation of 
the blade-element forces and moments in a rotating frame. 
The “Technion Model” uses aerodynamic harmonic func- 
tions as forcing terms to the equations of motion expressed 
in a non-rotating frame. The two models have been config- 
ured to represent the AH-64 Apache helicopter (Fig 1 ), and 
an enhanced simulation capability for the AH-64 has be- 
come available as a result of this work. This paper presents 
an overview of both modeling approaches and a compari- 
son of the simulated responses against available flight-test 
data. This comparison of methodology based on a com- 
mon aerodynamic and flight-test data-base permits a good 
opportunity to assess the tradeoffs between simplicity and 
accuracy. 

Modeling Methods 

This section describes the two methods used for modeling 
the rotor. The airframe modules for both BEM AP and the 
Technion Model are essentially the same since they both 
use look-up tables based on the same set of wind-tunnel 
data. 

Overview of the Approach 
Used by the U.S. Army 

The structure of the Sikorsky GEN HEL main-rotor mod- 
ule (Ref 5) was chosen as the starting point for the 
new AH-64 main-rotor. This choice was made because 
the GEN HEL main-rotor was the only readily available 
blade-element type module usable for real-time operation 
(Ref 6). The equations for the coupled flap-lag dynamics 
of each blade were derived symbolically with the aid of 
the symbolic manipulation program MACS YMA (Ref 7). 
A Newtonian approach was used for the derivation and no 
simplifying assumptions were made other than the use of 
rigid blades. The equations were completely expanded to 


allow a close look at the effects of various toms, and to 
allow relative magnitude analysis on the terms that are of- 
ten ignored to ensure that they are indeed negligible for the 
flight condition being considered. The complete equations 
were retained for this work. 

The UH-60 specific equations in the GEN HEL main-rotor 
module were replaced with the newly derived equations 
and the UH-60 specific data in the module replaced With 
the corresponding AH-64 data. The MDHC model of the 
Apache, known as FLYRT (Ref 4) and obtained under con- 
tract to the U.S. Army, was then restructured to allow the 
upgrading of the model by replacing the main-rotor mod- 
ule. A few other modules also had to be upgraded to sup- 
port the new blade-element rotor. Finally, the input/output 
structure of the rotor module was revamped to allow inter- 
face with FLYRT instead of GEN HEL. The U.S. model of 
the AH-64 being used here is therefore a restructured and 
updated FLYRT employing a blade-element type rotor in- 
stead of the model’s original map-type main-rotor module. 
The original FLYRT with the Rotor-Map main-rotor mod- 
ule was recently validated by comparison with available 
flight data (Ref 8). The same flight data will be used in the 
present report for ease of comparison. 

Formulation of Rotor Equations; As mentioned pre- 
viously, the equations for the coupled flap-lag motion and 
the rotor force and moment were derived with the aid 
of the symbolic manipulation program MACSYMA. Fig- 
ure 2 depicts the coordinate systems used for the Newto- 
nian derivation. Note that a flap-lag-pitch hinge sequence 
has been used which does not include the pitch-lag cou- 
pling that exists with the flap-pitch-lag hinge sequence of 
the Apache. It was decided that the improvement in model 
accuracy afforded by the inclusion of the pitch-lag cou- 
pling does not justify the significant increase in equation 
complexity that results from the inclusion. 

First, the hub inertial translational and rotational veloc- 
ity and acceleration vectors were derived (rotating-shaft 
frame) based on the rotational and translational velocity 
and acceleration vectors at the aircraft C.G., and the rel- 
ative location of the C.G. and the rotor hub. Then, the 
velocity and acceleration vectors at a given blade-element 
were derived (rotating-shaft frame), as the sum of the hub 
motion and the local flapping and lead-lag. 

The coupled flap-lag equations of motion were then de- 
rived as a balance of inertial, aerodynamic, gravitational, 
and restraint (flap and lead-lag spring and damping terms) 
moments about the lead-lag and the flapping hinges. The 
inertial terms of these equations were compared with a 
previous, Lagrangian derivation by Chen (Ref 9) and 
shown to be a perfect match. The rotor force and moment 
vectors were also derived and MACSYMA was used to 
generate FORTRAN code corresponding to the new equa- 
tions. The aerodynamic forces for each blade-element are 
calculated using swept wing approximations (Ref 5), with 
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table look-up of lift and drag coefficients as a function of 
the local angle of attack and Mach number. The look-up 
tables were constructed based on data available in the “Air 
Vehicle Technical Description Data for the AH-64A Ad- 
vanced Attack Helicopter” (Ref 10). Since a good model 
of the lead-lag dampers was not available, flapping and 
lead-lag spring and damping terms were included in the 
equations as a temporary alternative. The values of these 
parameters were provided by McDonnell Douglas Heli- 
copter Company. 


The rotor force and moment vectors were calculated by 
numerically summing the elemental forces and moments, 
first over each blade and then over all the blades. Spe- 
cial attention was given to the transfer of moments through 
each hinge as components along the two axes not aligned 
with the axis of the hinge. These moments are often ig- 
nored by assuming that hinges do not transfer any- mo- 
ments, which, of course, is only true of the components of 
moments aligned with the axes of the hinges. This may be 
seen if we look at the detailed derivation of the rotor mo- 
ments from a blade-element down to the hub. Let A Ft, 
AFr, and A Fp be the tangential, radial, and perpendic- 
ular components of the forces on a typical blade-element 
(Fig 2). Then, the elemental force in frame 1 may be writ- 
ten as 


AF/i = 

cos 8 
— sin 8 

sin 8 
cos 8 

0 ‘ 
0 
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f -AFr 1 
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0 

0 
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l -AFr J 


-A F t cos 8 + A Fr sin 6 
A Ft sin 6 + AFr cos 8 
-A F P 


( 1 ) 


The moment arm from the outer hinge (lead-lag in our 
case) to the blade-element may be written, in frame 1, as 



Therefore, the elemental moment at the outer hinge (lead- 
lag), which is the cross product of the moment arm and the 
elemental force, is 
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The elemental moment at the inner hinge (flapping) is the 
sum of the moment transferred through the outer hinge 
(lead-lag) and the moment due to the shear face at the 
outer hinge. Therefore, the elemental moment at the inner 
hinge, in frame 2, may be written as 
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Finally, the elemental hub moment, in the rotating-shaft 
frame, may be written as the sum of the moment trans- 
ferred through the inner hinge and the moment due to shear 
face at the inner hinge. Therefore, the elemental moment 
at the hub, in the rotating-shaft frame, may be written as 

A M hub, r» 
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Only the second vector of the right-hand side, which rep- 
resents the moment of the shear face at the inner hinge 
multiplied by the inner hinge offset, is usually considered 
in derivations (Ref 5). However, experimentation with in- 
troducing the appropriate extra terms in GEN HEL (for the 
UH-60 lag-flap-pitch sequence) (report by M. H. Mansur 
to be published as an Army-NASA Technical Memoran- 
dum) has shown that the extra terms may also be signifi- 
cant They were therefore retained for this work. 


Overview of Approach Used at the Technion 


Researchers of the Faculty of Aerospace Engineering at 
the Technion have been developing rotorcraft flight me- 
chanics simulation models for the last 10 years (Refs 11- 
15). The initial approach to rotor modeling was based on 
a Tip-Path Plane approach, or more accurately, taking into 
account only the constant and first harmonic of the forcing 
terms in the rotor equations of motion and the response 
variables. Recently, the model has been extended and 
is capable of taking into account higher harmonics also. 
More details appear in the subsection on the formulation 
of the rotor equations. 

While the rotor represents the most important element of 
the helicopter, a balanced model requires an appropriate 
(accurate enough) description of the contribution of other 
elements of the vehicle. The method of dealing with the 
contributions of these elements in trim calculations was 
first described in Ref 14. This method was lata extended 
to include maneuvers and stability calculations as well. A 
very brief description of the method of calculating these 
contributions, which include fuselage, tail-rotor, etc., will 
be presented for the sake of completeness. 

1. Fuselage: Inertia contributions are dealt with in an 
accurate manna and look-up tables, based on wind 
tunnel tests, are used to calculate the aerodynamic 
contributions. 
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2. Tail Rotor The dynamics include only flapping and 
the calculations a re similar to the main rotor. First- 
order interference effects between the vertical fin and 
the tail rotor are included. 

3. Vertical fin, wings, stabilizer. The calculations of the 
contributions of the aerodynamic surfaces are based 
on look-up tables. Corrections for side-slip and as- 
pect ratio are included. 

4. External stores (may include rockets, missiles, fuel 
tanks, etc.): The inertia contributions are added to 
the fuselage. Look-up tables are used for the aerody- 
namic calculations. 

Formulation of Rotor Equations; The Technion rotor 
model is an extension of the model that was described in 
Refs 11-13. Only a brief description is presented here. 
More details will appear in forthcoming publications. 

The blade’s equations of motion are derived using La- 
grange’s equation. The general form of these equations 
is 

U]{a} + Q[B]{d} + 0 2 (C]{<r} 

= {Q/}+{OJ (6) 

where {ff} is the vector of unknowns: 

{a} r = </J,Git,0e> (7) 

and p, 0 R , and Q c define the angle of flapping, elastic pitch 
at the blade root, and elastic pitch variations along the 
blade respectively. 

[A], [B], and [ C] in Eq (6) are square matrices of or- 
der 3 that include all the design details of the rotor and 
the blades. {Q/} is the forcing vector that includes all the 
effects except aerodynamics. {Q^} is the forcing aerody- 
namic vector, defined as 

{QaY = < QfiA, QtRA, QttA > ( 8 ) 

and QfiA, Q»ra . and Q» k a are generalized aerodynamic 
loads. The expressions for these loads are obtained by ap- 
plying the principle of virtual work. A special ordering 
scheme is applied to simplify the equations. 

At each time step during the simulation the generalized 
aerodynamic loads are calculated at a finite number of az- 
imuthal locations. Then, by using a Fast Fourier Trans- 
form procedure the generalized aerodynamic loads are ex- 
pressed in the following harmonic form: 


Ni 

Q$a = Mpo + TJMm cos(;V 0 
+M^s, sin(;V>)] 

Qbra - cos (jip) 

1 

+ MeRs t sin(;'i/>)] (9) 

N i 

QttA = Affeo + y '[MetCj cosQVO 
>" i 

+M itS , sin (jip)] 

where ip is the blade’s azimuth angle and Nj is an input to 
the computer program. As Nj increases the accuracy in- 
creases, but, at the same time, the required computer time 
is increased as well. 

After the three coupled equations of motion of a single 
blade are obtained, a multiblade coordinate transformation 
is applied. As a result, each of the three variables p, 0 R , 
and 8c is replaced by four new unknowns associated with 
the multiblade coordinates for a four-bladed rotor. ()o, 
Oc.Os, and ( Xy/ 2 . Thus, Eq (6) for each blade is replaced 
by the following multiblade equation: 

+[D m ]{pq}+[E m ]{PC) 

+(F m ]{PC}+{f m } 

= M a (10) 

where r is the vector of rotor variables: 

{t) t = < Po,0ro,O<*>,Pc,0rc,QcC,Ps, 

GrS, ®eS, Pn/2 , OrN/2 , O'N/2 > (H) 

and {pg} is the vector of angular rates of the hub and 
{PC} is the vector of collective and cyclic pitch inputs 
to the main rotor. 

{Ma} in Eq (10) is the vector of aerodynamic loads after 
transformation to the multiblade coordinates. This vec- 
tor is a function of the harmonic coefficients of Eq (9). It 
should be pointed out that the aerodynamic calculations 
include dynamic-inflow effects (Ref 16) and aerodynamic 
interference between the rotor and the fuselage (Ref 17). 

Figure 3 represents the aerodynamic flapping moment at 
the blade root in trimmed horizontal flight at an airspeed 
of 100 knots. Two cases, one for Nj = 1 and the other for 
V ; - = 4 , are compared. In the case of N } - = 1 , only the 
first harmonic terms are considered in the blade flapping 
motion. The figure shows (Nj = 1 actual) that even with 
only the first harmonic terms included in the blade dynam- 
ics, the actual flapping moment includes higher harmonic 
components. The first harmonic approximation to this ac- 
tual flapping moment is also shown in the figure {Nj = 1 
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approximate). It may be clearly seen that fairly large devi- 
ations exist between the first harmonic approximation and 
the actual flapping moment 

In the case Nj - 4 , the first four harmonics are included 
in the blade dynamics. As may be seen, the actual flap- 
ping moment again includes high harmonic components, 
even above the first four harmonics. The presence of har- 
monics above the first four is indicated by the slight de- 
viation between the actual flapping moment and its first- 
four harmonics approximation (Nj = 4 actual compared 
to Nj = 4 approximate). The close match between the two 
curves, however, suggests that including harmonics higher 
than the first four will result in only negligible changes. 
The figure also shows that the inclusion of higher harmon- 
ics in the blade dynamics results in differences in the ac- 
tual flapping moment, as seen when the actual flipping 
moment for the case of including the first four harmonics 
in the blade motion (Nj = 4 ) is compared with the case 
of including only the first harmonic (Nj =1). 

An important and interesting result to note is that despite 
the difference in the actual flapping moment between the 
cases Nj = 1 and Nj = 4 , the first harmonic approxima- 
tion to the actual flapping moment for the case Nj = 4 
is almost identical to the first harmonic approximation to 
the actual flapping moment for the case of Nj - 1 . This 
explains why the first harmonic approximation gives such 
good results in flight mechanics problems where frequen- 
cies of only up to 1/rev are of interest. This fact is further 
illustrated in Fig 4 which depicts the total lift transferred 
from a blade to the hub. 

In order to increase the model efficiency in trim calcu- 
lations, the direct integration with respect to time is re- 
placed by a solution of a fairly large nonlinear periodic 
problem. This problem is solved using a new method of 
obtaining numerical solutions for highly nonlinear peri- 
odic problems (Ref 18). 

Thus, the present Technion Model offers a very conve- 
nient way of changing the accuracy of the rotor model. 
By choosing values of Nj between 1 and very large, one 
is capable of “moving continuosly” between a “tip path 
plane” approach and a classical blade-element straightfor- 
ward integration with respect to the azimuth. 

Results 

The ability of each model to correctly simulate the flight 
mechanics of the AH-64 was determined by comparisons 
with available flight test data. The comparisons were per- 
formed by trimming the simulation models to the flight 
condition being considered and driving the models with 
the recorded flight test control inputs in all axes. To pre- 
vent the trim discrepancies between the models and the 
actual flight vehicle from affecting the dynamic response 
comparisons, calculated variations from the flight-test- 


control trim values were used in constructing the simula- 
tion inputs. In other words, the simulation control inputs 
consisted of the sum of the simulation trim values and the 
flight test control variations from trim. Both the on- and 
the off-axes responses of the models were looked at and 
compared with flight-test data. This paper, however, con- 
centrates on the on-axis responses for die sake of brevity, 
and includes only one off-axes case for completeness. The 
difficulties of matching coupled responses are also briefly 
discussed. 

Data obtained under contract from MDHC were used for 
the model comparison work presented here. In addition to 
being the best data available, the use of the MDHC data 
allows a comparison of the responses of the new mod- 
els with the MDHC model FLYRT, since the same data 
were used for validating FLYRT as outlined in Ref 8. 
The reader is referred to that document for information 
regarding flight conditions, processing, and consistency 
checks/corrections performed by MDHC. 

Even though the Technion Model and BEMAP are com- 
pletely different in origin and design, they are based on the 
same set of wind-tunnel and flight-test data. Furthermore, 
they both employ the same basic modeling approaches in 
all the modules except the main rotor. The major differ- 
ences between the two rotor-modeling approaches docu- 
mented above may be summarized as follows: 

1 . BEMAP uses a direct blade-element approach which 
considers all the harmonics, whereas the Technion 
Model includes only the first harmonic since the first 
harmonic was shown to almost solely dominate the 
aircraft response for flight mechanics applications. 

2. The Technion Model uses the Pitt/Peters 
(Ref 19) dynamic inflow model while BEMAP can 
use either Pitt/Peters or the “extended” Howlett in- 
flow model described by Ballin (Ref 20). The two 
inflow models are quite similar, as confirmed by com- 
paring BEMAP responses using Pitl/Peters with BE- 
MAP responses obtained using the extended Howlett 
model. 

3. The airfoil tables used for both models are the same 
(Ref 10). However, whereas BEMAP incorporates 
look-up tables for calculating section forces, the 
Technion Model uses polynomial fits of the data in 
the rotor derivations. 

4. The Technion Model does not consider the lead/lag 
degree of freedom whereas BEMAP does. This, 
however, does not necessarily help BEMAP because 
of the unavailability of a good model of the elas- 
tomeric lead-lag dampers used on the Apache. 

Static Validation 

Static validation refers to the comparison of the equilib- 
rium trim conditions of the aircraft and simulation models. 
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Figures Sa-f show that both models simulate the trim con- 
trols and the attitudes of the helicopter quite well across 
a range of airspeeds. The trim comparisons do not point 
to a clear choice as far as modeling the rotor is concerned 
since the models alternately come closer to the flight data 
as is seen from the figures. For example, figure Sa shows 
that BEMAP consistently underestimates the trim pitch at- 
titude of the aircraft, up to a maximum error of about 2 de- 
grees, above an airspeed of around 40 knots. Also, as fig- 
ure Sb indicates, the direct blade-element model overesti- 
mates the collective input required to trim by slightly over 
half an inch in the range of 40 to 80 knots whereas the har- 
monic model does quite well up to around 120 knots. In 
calculating longitudinal cyclic required to trim, however, 
BEMAP duplicates the trends in the flight data quite accu- 
rately whereas the Technion Model is slightly off at lower 
speeds, as may be seen in figure 5c. The same is true for 
the lateral cyclic required to trim, as may be seen in fig- 
ure 5e. Finally, figures 5d and 5f show that the two models 
are quite comparable as far as calculating the trim roll an- 
gle and trim directional control are concerned. Therefore, 
until further work has detrmined the cause of the discrep- 
ancies in each model, neither can be judged better or worse 
than the other. 

Dynamic Validation 

Dynamic validation refers to the comparison of the dy- 
namic response characteristics of the aircraft and simula- 
tion model following a control input away from trim. Dy- 
namic response is much more difficult to simulate accu- 
rately than trim because the random conditions present at 
the time of the flight tests, such as wind and turbulence, af- 
fect the responses of the aircraft but are difficult to record 
or model. Post processing schemes to remove the random, 
uncorrelated responses may be used to reduce the severity 
of the problem (Ref 21), but none woe attempted here. 

The on-axis responses of both models are good (only 
DASE-off data were considered for this work). Con- 
centrating on the slopes of the various coplotted lines in 
Figs 6-11, rather than the absolute values, we see that in 
almost all the cases both models simulate the overall re- 
sponse fairly accurately. The mismatches usually seem 
to be in the form of a bias starting from a mismatch of 
trim conditions. Both models seem to overpredict the rate 
of onset of the accelerations by about the same amount 
This is interesting given the fact that the rotor is mostly re- 
sponsible for the accelerations and that the two models use 
completely different rotor modules. The secondary effects 
not modeled by either model, as opposed to the differences 
between the two models, are most likely responsible for 
the variations from the actual aircraft response. These sec- 
ondary effects include compressibilty, reverse flow, and 
tip losses. 

It is also interesting to note that die Technion Model seems 
to correlate better in the lateral axis, whereas BEMAP 
seems to correlate better in the longitudinal, regardless 


of flight condition. The lateral responses of the BEMAP 
model at both hover and 80 kts (Figs 6 and 7) is seen to be 
less damped than the Technion Model whereas the longitu- 
dinal responses of the BEMAP model (Figs 8 and 9) seem 
to be overall closer to the flight data than the Technion 
Model. One possible reason for this may be the lead-lag 
degree of freedom. The Technion Model does not consider 
the lead-lag and it appears that, at least as far as lateral re- 
sponse is concerned, not considering the lead-lag may be 
better than considering it without a representative model 
of the dampers. 

Figures 10 and 11 depict the directional responses of the 
models to directional doublets at hover and 80 kts, respec- 
tively. As may be seen, in both cases BEMAP trims with 
excessive right pedal and fails to correlate with flight data 
as well as it does in the other axes. The Technion Model, 
on the other band, does quite well at hover and somewhat 
better than BEMAP at 80 kts. Introduction of a good lead- 
lag damper may improve BEMAP in this axis as well. 

The off-axis responses of both models are generally less 
accurate than the on-axis. Figure 11 shows the lateral re- 
sponse to a longitudinal doublet at hover. As may be seen, 
both models trim quite well but do not seem to duplicate 
the dynamic response very accurately. This is to be ex- 
pected since the off-axis response is dominated by cou- 
pling and secondary effects which are difficult to model 
accurately. 

Concluding Remarks 

Both the BEMAP and the Technion Models are still in the 
development phase and are being continually updated to 
improve their accuracy and sophistication as necessary. 
Secondary effects such as compressibility, reverse flow, 
and tip loss have not yet been sufficiently incorporated in 
either model. As seen from the results presented, however, 
even at this stage both models are capable of simulating 
the response of the AH-64 helicopter quite closely. Cur- 
rently, flight tests are being conducted at the Army Avia- 
tion Flight Activity (AEFA) to provide additional data to 
resolve the validation discrepancies. 
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Fig 1: AH-64 Apache Advanced Attack Helicopter. 
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Fig 2: Coordinate system definitions for the Blade-Element Rotor derivations. 
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Fig 3: Aerodynamic flap moment versus azimuth angle. 



Fig 4: Total lift force versus azimuth angle. 
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Fig 6: On-axis response to lateral doublet left, at hover. Fig 7 : On-axis response to lateral doublet right, at 80 knots. 
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